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ABSTRACT 

We present Herschel PACS 100 and 160 pm observations of the solar-type stars a Men, HD 88230 and HD 210277, which form part of the FGK 
stars sample of the Herschel Open Time Key Programme (OTKP) DUNES (DU st around A/Earby S tars). Our observations show small infrared 
excesses at 160 pm for all three stars. HD 210277 also shows a small excess at 100 pm, while the 100 pm fluxes of a Men and HD 88230 agree 
with the stellar photospheric predictions. We attribute these infrared excesses to a new class of cold, faint debris discs, a Men and HD 88230 
are spatially resolved in the PACS 160 pm images, while HD 210277 is point-like at that wavelength. The projected linear sizes of the extended 
emission lie in the range from ~ 115 to < 250 AU. The estimated black body temperatures from the 100 and 160 pm fluxes are < 22 K, while 
the fractional luminosity of the cold dust is L dust IL+ ~ 10 -6 , close to the luminosity of the Solar-System’s Kuiper belt. These debris discs are the 
coldest and faintest discs discovered so far around mature stars and cannot easily be explained by invoking ’’classical” debris disc models. 

Key words. - Stars: planetary systems: planetary discs -Stars: planetary systems: formation - Stars: individual: a Men (HIP 29271) - Stars: 
individual: HD 88230 (HIP 49908) - Stars: individual: HD 210277 (HIP 109378) 


1. Introduction 

Debris discs are tenuous structures associated with main se- 
quence stars formed by second generation dust, which has re- 
sulted from the collisions of solid bodies continuously supplying 
the circumstellar environment with small dust particles. This for- 
mation sequence is inferred from the lifetime of the dust grains 
against destructive collisions, Poynting-Robertson drag and ra- 
diation pressure, which is much shorter than the ages of the host 
stars. These discs are visible in reflected light at optical wave- 
lengths and in thermal radiation at mid-/far-IR and submillimeter 
wavelengths (Aumann et al., 1984; Backman & Paresce, 1993). 
General disc characteristics are grain black body temperatures of 
~ 50 - 100 K, fractional luminosities / > ~ 10 -5 , and 

radii from less than 10 AU to several times 100 (e.g. Absil et al. 
2006, Su et al. 2005, Trilling et al. 2008). Debris discs are con- 
sidered analogues of the Solar System asteroid and Kuiper belts, 
although their luminosities are usually more than 100 times the 
Kuiper belt level of L dust /L Su n ~ 10 -7 - 10 -6 (Stern 1996, 
Vitense et al. 2010). The sensitivity of the 3.5 m Herschel far- 
infrared space telescope (Pilbratt et al. 2010) with its instrument 
PACS (Poglitsch et al. 2010) offers the possibility of characteris- 
ing colder (~ 30 K) and fainter (. L dus t /L* few times 10 -7 ) debris 
discs with spatial resolution ~ 60 AU (FWHM) at 10 pc, i.e., 
true extra-solar Kuiper belts. 

DUNES is a Herschel OTKP designed to detect and char- 
acterise extra- solar analogues to the Kuiper belt around main 
sequence FGK nearby stars (Eiroa et al. 2010). In this letter 
we present the results for three stars from the DUNES sam- 
ple: a Men, HD 88230, and HD 210277 as clear examples of 

* Herschel is an ESA space observatory with science instruments pro- 
vided by European-led Principal Investigator consortia and with impor- 
tant participation from NASA. 


Table 1 . Stellar properties. 


Star 

a Men 

HD 88230 

HD 210277 

Sp. Type 

G5V-G7V 

K6V-M0V 

GOV, G7V-G9V 

Distance (pc) 

10.2 

4.9 

21.6 

L(L 0 ) 

0.85 

0.15 

1.10 

T«//(K) 

5590 

3850 

5540 

Age (Gyrs) 

5.5 

6.6 

6.9 


the advantages offered by Herschel observations: they trace a 
new class of cold, T < 22 K, spatially resolved debris discs with 
very low fractional luminosities. These discs have remained un- 
observed by previous far-IR and submillimeter studies. Table 1 
gives some properties of the stars. Ages are based on the log 
R' HK activity index and have an uncertainty of 60% (Mamajek 
& Hillenbrand 2008). HD 210277 hosts a Jupiter-like planet 
(Marcy et al. 1999). This star and a Men have faint stellar com- 
panions, but neither the measurements nor the photospheric pre- 
dictions are affected by them. Eiroa et al. in prep will present a 
full discussion of the stars and the general results of the DUNES 
survey. 

2. Observations and data reduction 

a Men, HD 88230 and HD 210277 were observed with PACS 
100/160 pm in scan map mode. For each star, two scans at po- 
sition angles 70° and 110° were carried out, each scan consist- 
ing of 10 legs with separation of 4", length of 3' and medium 
speed of 20"/s. Table 2 gives the scan identification numbers 
(Obs. ID) and the total duration of the observations (OT). a 
Men was observed twice in order to increase the S/N ratio. Data 
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Table 2. Log of the PACS 100 pm and 160 p m observations. 


Star 

HIP 

Obs. ID 

OT (sec) 

a Men 

29271 

1342203297/8 

1116 

a Men 


1342216043/4 

2244 

HD 88230 

49908 

1342210610/1 

4500 

HD 210277 

109378 

1342211126/7 

4500 


reduction was made using the Herschel Interactive Processing 
Environment (HIPE) version 7.2. The individual scans were mo- 
saiced to produce the final image at each band. To check the con- 
sistency of the reduction and analysis (particularly regarding the 
effect of correlated noise), mosaics were produced at both the na- 
tive 372 for 100 /mi (green) and 674 for 160 p m (red), as well as 
super-sampled pixel scales 1.0" (green) and 2.0" (red), the latter 
being the default pixel size in HIPE for this type of data. The 
images with the native pixel size scales avoid, at least partly, the 
correlated noise in the PACS image^](see also Fruchter & Hook, 
2002). A high-pass filter was used to remove large scale back- 
ground emission from the images, with filter widths of 15" and 
25" in the green and red channels. To prevent the removal of any 
extended structure near the stars, regions where the sky bright- 
ness exceeded a threshold value in the image defined by the stan- 
dard deviation of all the positive pixels ( S > 10 ~ 6 Jy/ pixel) 

were masked from that process. Absolute flux calibration un- 
certainties are ~3% and ~5% for the green and red bands (see 
technical note below). 

3. Results 

Table [3] gives the J2000.0 optical equatorial coordinates of the 
stars as well as their 100 pm peak positions corrected for the 
proper motions of the stars (van Leeuwen, 2007). Offsets be- 
tween the optical and the PACS positions (column 4 of Table [3]) 
are within ~ 1.5<x the Herschel pointing accuracy of 274 in this 
observing mode (Sanchez-Portal, private communication). 

Fig. 1 shows the 100 and 160 pm images and isocontour plots 
of the stars, a Men and HD 88230 are point-like at 100 pm, 
FWHM ~ 673 x 6 75, while both stars are resolved at 160 /mi, 
with angular sizes/position angles ~ 18"/52° and ~ 2375/45°, re- 
spectively. The extended emission appears clearly asymmetric in 
HD 88230, with the star located at the North-Eastern side. None 
of the objects are resolved in the direction orthogonal to the ex- 
tended emission. HD 210277 is unresolved at both wavelengths. 

PACS fluxes (Table [3]) have been estimated using circu- 
lar and rectangular aperture photometry, taking special care to 
choose the reference background region due to the presence of 
field objects. Specifically, measurements of HD 210277 take into 
account the presence of the bright object located at ~16" North- 
East from the star; in the case of this star we have, in addition, 
carried out PSF photometry using the DAOPHOT software pack- 
age. The PSF photometry fluxes are F(100) = 8.4+0.3 mJy and 
F(160) = 14.3+0.4 mJy (errors are those of the PSF fits), which 
are consistent with the aperture photometry estimates. Errors 
have been estimated using a variety of methods, including cir- 
cular and rectangular boxes at different positions in the nearby 
fields. Sky sizes for the errors estimates are equal to the area of 
the aperture used for the photometry for all three stars and both 
bands; in particular the sizes of the extended emission around a 
Men and HD 88230. 

PACS fluxes have been compared to predicted stellar photo- 
spheric fluxes (Table [3} using Gaia/PHOENIX models (Brott & 

1 Technical Note PICC-ME-TN-037 in http://herschel.esac.esa.int 


Hauschildt, 2005), with the stellar parameters as given in Eiroa 
et al. (in prep). Fig. 1 shows the spectral energy distributions 
(SEDs) of the stars, where PACS fluxes are plotted together with 
optical, near-IR, IRAS, AKARI, and Spitzer / MIPS and IRS data; 
in addition, the best ^-square photospheric fit is shown. To assess 
the presence of an excess at 100 and/or at 160 pm we require that 
the observed fluxes, Epacs, exceed by at least 3<x the predicted 
photospheric fluxes, F * (xx = (Fpacs~ F*)/o- A ). No excesses are 
detected at 100 pm for a Men and HD 88230, while it is seen in 
HD 210277. All three stars do show 160 pm excesses. We also 
note that the SED slopes from 100 to 160 pm are a = - 0.5+0.7 (a 
Men), a = - 0.7+0.6 (HD 88230), and a = 0.8+0.8 (HD 210277), 
which clearly differ from the expected Rayleigh- Jeans behaviour 
(a = - 2.0) of a stellar photosphere in this wavelength regime. 

3.1. Comments on the PACS images 

All three fields show ~ 2-3 red sources per square arcmin, which 
are likely to be background galaxies. In particular, there is a faint 
peak at ~13" towards the SW from HD 88230 and a bright one at 
~16" towards the NE from HD 210277. We have consulted the 
NASA/IPAC extragalactic database to search for counterparts 
without finding any association. Although we cannot firmly ex- 
clude a coincidental alignment or contamination of a background 
source(s) in the line of sight of our stars, we think that it is un- 
likely for the Herschel sources presented here due to the close 
correlation between the optical and Herschel positions, and the 
photospheric predictions and the estimated Herschel 100 pm 
fluxes. In fact, following in a first approach the source counts by 
Berta et al. (2010), the average density of extragalactic sources 
with fluxes ~ 6-7 mJy and ~ 12 mJy - i.e., the measured excesses 
at 160 pm from a Men/HD 88230 and HD 210277, respectively 
(Table 3) - are 2/arcmin 2 and 0.7/arcmin 2 . Thus, given the opti- 
cal/160 pm offsets (Fig. 1) the a priori probability of an acciden- 
tal alignment is very small, clearly smaller than 5%. We refer to 
a future paper (del Burgo et al., in prep.) for a detailed study on 
source contamination within the DUNES fields. 


4. Analysis 

The small PACS excesses above the photospheric fluxes are in- 
terpreted as due to cold debris discs around the stars. Black body 
dust temperatures, 7d US t, can be estimated from the 100 and 160 
pm excess fluxes for HD 210277; in the case of a Men and HD 
88230, an upper limit for 7d us t can be calculated taking the 100 
pm flux as 3<x statistical noise. 7d US t for the three stars is < 22 
K (Table 4); the corresponding inner radii of discs with black 
body grains at this T^ ust , considering the luminosity of the stars, 
are also given in Table 4. Fig. 1 shows the excellent agreement 
between the combined SED of the stellar photospheric fits plus 
22 K black bodies and the observed SEDs. 

Deconvolution of the images can be used to estimate the true 
size of the resolved discs in a Men and HD 88230. Our method 
first removes the photospheric contributions from each image 
by subtracting a PSF with a peak scaled to the predicted photo- 
spheric flux level. The PSF model uses a Bootis images, rotated 
to match the roll angle of the telescope during the observations of 
the DUNES stars. After star subtraction, the images are decon- 
volved using both modified Wiener and Richardson-Lucy algo- 
rithms. The noise model takes into account that the main contrib- 
utor is the telescope thermal emission. Both algorithms produce 
consistent results although with different noise patterns. Fig. 2 
shows the 160 pm star- subtracted and the Wiener-deconvolved 
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Table 3. Optical positions of the stars and of their PACS 100 identified counterparts. Observed PACS fluxes with l<x statistical 
errors (FpacsX and predicted photospheric fluxes (E*). Flux units are mJy. 


Star 

Optical position 
(J2000.0) 

PACS 100 yum position 
(J2000.0) 

Offset 

PACS 100 yum 
Epacs E* xioo 

PACS 160 ji m 
EpACS E* X\60 

a Men 

06 10 14.47 -74 45 11.0 

06 10 14.53 -74 45 11.1 

0.3 

17.8+1.3 

18.0 

-0.2 

14.4+2.0 

7.0 

3.7 

HD 88230 

10 11 22.14 +49 27 15.3 

10 11 21.88 +49 27 17.2 

3.2 

22.5+0.9 

25.7 

-3.6 

16.0+1.7 

10.0 

3.5 

HD 210277 

22 09 29.87 -07 32 55.2 

22 09 30.08 -07 32 54.1 

3.3 

8.5+1.0 

4.6 

3.9 

12.4+1.6 

1.8 

6.6 



Fig. 1 . PACS images, isocontours and SEDs of a Men/HIP 29271 (left), HD 88230/HIP 49908 (middle) and HD 210277/HIP 109378 
(right). Images: 100 yum (left), 160 yu m (right). North is up and East to the left. Isocontours: 100 yum contours are in green colour 
while the 160 yum ones are in red. a Men: 100 yu m contours (10, 20, 40, 80 and 90 % of the flux peak); 160 yum contours (20,40,60,80, 
90 % of the flux peak) HD 882308: 100 yum contours ( 10, 20,30,40,60, 80, 90 % of the flux peak); 160 yum contours (20, 30, 40, 
50, 60, 80, 90 % of the flux peak) HD 210277: 100 yum contours ( (30, 40, 50, 60, 70, 80, 90 % of the flux peak); 160 yum contours 
(50, 60, 70, 80, 90 % of the flux peak). The lowest contour in all cases is ~ 3 cr. The optical position of the stars are indicated by the 
symbol in the isocontour plots; a segment indicates the projected linear sizes at the distance of each star. SED plots: black line 
is the photospheric fit while the blue line is the photosphere plus a 22 K black body. 


images of a Men and HD 88230; estimated sizes of the decon- 
volved sources are 16" and 21", respectively. Table 4 gives the 
linear sizes of the semi-major axes (~3<x contours) from both the 
original and deconvolved PACS 160 yum images. In the case of 
HD 210277, the value in Table 4 corresponds to an upper limit 
of 12", i.e., approximately the 160 /rm beam size. The compar- 
ison of both the directly observed and deconvolved semi-major 
axes to the estimated radii from T^ust clearly indicates that the 
observed discs (assuming as a first approach symmetric discs) 
are smaller than the expected sizes for black body discs. 

Dust fractional luminosities, / (Table 4), can be obtained by 
taking the 160 jim fluxes and assuming the star temperatures 
given in Table 1 and 7d us t = 22 K (Beichman et al., 2006). The 


estimated / values are of the order of 10 -6 , close to the Kuiper 
belt level. Essentially similar / values are obtained taking the 
maximum wavelength corresponding to 22 K, i.e., 230 /im, and 
its expected flux extrapolating from the one measured at 160 yum. 


5. Discussion 

The debris discs in this work are a new class of discs charac- 
terised by an excess at 160 /rm, little to no excess at 100 yum, 
and no excess emission at shorter wavelengths. These discs are 
the coldest and least luminous ones known to date; they are sig- 
nificantly colder and fainter than other observed DUNES discs 


4 


C. Eiroa et al.: Herschel discovery of a new class of cold, faint debris discs 




Fig. 2. 160 n m star- subtracted (left) and deconvolved (right) im- 
ages. Up: cr Men/HIP 29271. Contours: 10%, 20%, 40%, 80% of 
the peak. Bottom: HD 88230/HIP 49908. Contours: 20%, 40%, 
80% of the peak. North is up and East to the left. Field size is 
60" x 60". 

Table 4. Black body dust temperatures from the far-IR excesses, 
r dust and the corresponding estimated radii. Observed 3 cr linear 
sizes (semi-major axes) of the extended emission at 160 p m in 
the original and deconvolved images. / is the dust fractional lu- 
minosity. 


Star 

Tdus, 

(K) 

Size 

Est. 

(AU) 

Size 

Orig. 

(AU) 

Size 

Deconv. 

(AU) 

/ 

a Men 

<22 

>147 

92 

81 

9.7xl0 _v 

HD 88230 

<22 

>62 

56 

51 

1.6xl0“ 6 

HD 210277 

22 

160 

<130 


5.4xl0 -6 


(e.g. Liseau et al. 2010, Marshall et al. 2011). Their interpreta- 
tion poses significant challenges. 

The shapes of the SEDs suggest that the dust is located in 
a ring with a larger inner void. We have seen that the observed 
disc radii are smaller than those implied from T dust (Table 4). In 
other words, T dust is smaller than the black body temperature, 
Tbb ~ 25-30 K, given the observed sizes and stellar luminosi- 
ties. Thus, we need to reconcile these seemingly contradictory 
results. Since one would reasonably expect a broad grain size 
distribution, which would include small grains hotter than those 
indicated by the estimated T dust , the question is how to make the 
dust cold, while allowing a size distribution, but still keeping the 
disc radii within the observed ones. 

The low temperatures of grains require them to be large and 
highly reflective. The required albedo can be derived from r dust 
and Tbb- Taking T dust = 22 K as representative for the three ob- 
served discs, and r bb = 25-30 K, the albedo would be >50%. For 
comparison, surfaces of trans-Neptunian objects contain signif- 
icant amounts of ice (Barucci et al. 2011) and many, especially 
large ones (e.g. Pluto, Eris), have albedos in excess of 50% (e.g. 
Vitense et al. 2010). Accordingly, Pluto has a surface brightness 
temperature ~ 10 K below the black body value (Gurwell et al. 
2010). It is natural to expect that dust released from the surfaces 
of such objects would have similar properties. 


In addition, one has to explain why small grains in the cold 
discs are depleted. One possibility is to assume a very low dy- 
namical excitation of dust-producing planetesimals, so that discs 
would be devoid of small particles (Thebault & Wu, 2008). The 
reason is that low collision velocities between large grains, un- 
affected by radiation pressure, create an imbalance between the 
rates at which small grains are produced (low) and destroyed 
(high). Low collision velocities are compatible with low orbital 
velocities at the large radii of the cold discs. Besides, since the 
surface density of the solids far from the central stars is also 
low, planetesimal accretion scenarios predict very long accre- 
tion timescales (Kenyon & Bromley, 2008). Thus, large plan- 
etesimals that would excite the discs may have failed to grow. In 
fact, low dynamical excitation has been inferred for other large 
debris discs, e.g. HD 207129 (Lohne et al 2011). 

The tenet that the cold discs are probably in a low dynam- 
ical excitation state would be difficult to reconcile with the ex- 
istence of large planets in the discs, since they would stir the 
discs too strongly. The emission around a Men and HD 88230 is 
asymmetric, which might be suggesting the presence of a giant 
planet as in Fomalhaut (Kalas et al.2008). The problem can be 
mitigated if the planets are in nearly-circular orbits or the plan- 
etesimals have low eccentricities, as suggested for Fomalhaut 
(Chiang et al. 2009). HD 210277 hosts a planet at 1.1 AU with 
e p = 0.47; however, based on the formulae by Mustill & Wyatt 
(2009), a stirring front from such a planet should not be able to 
reach the -130 AU-sized disc on Gyr timescales. 

Another point is the origin of the large inner voids in the cold 
discs. These could either be due to clearing by planets or may 
reflect the accretional and collisional history of primordial discs. 
In the latter case, the observed radii are those at which solids 
could reach “right” sizes, “right” degree of dynamical excitation, 
and/or were able to survive over the stellar age, to produce the 
observed emission. A detailed analysis of the nature of the cold 
discs is beyond the scope of this paper, and we defer to Krivov et 
al. (in prep.) where possible scenarios will be discussed in detail. 

6. Conclusions 

We have presented Herschel PACS observations of three stars of 
the OTKP DUNES sample. The observations reveal a new class 
of debris discs with fractional luminosities close to the Solar- 
System Kuiper’s belt, but are colder and larger. These discs are 
a challenge to current models explaining debris discs around 
mature solar-type stars such as either the usual collisional- 
dominated disc scenario or low dynamical excitation discs. 
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